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Semiconductor Bloch equations, in their extension including the spin degree of freedom of the 
carriers, are capable to describe spin dynamics on a microscopic level. In the presence of free holes, 
electron spins can flip simultaneously with hole spins due to electron-hole exchange interaction. 
This mechanism named after Bir, Aronov and Pikus, is described here by using the extended semi- 
conductor Bloch equations [l| and considering carrier-carrier interaction beyond the Hartree-Fock 
truncation. As a result we derive microscopic expressions for spin-relaxation and spin-dephasing 
rates. 

PACS numbers: 67.57.Lm,72.25.Rb,42.50.Md 
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Semiconductor Bloch equations (SBE) are a well es- 
tablished concept to describe the dynamics of carriers in 
a semiconductor or quantum structure by a scalar light 
field. 0, S H It has been used successfully in modelling 
the time evolution due to carrier-carrier interaction on 
different time scales including the coherent and the re- 
laxation regime. Thus, SBE have become the domi- 
nating tool in the theory of semiconductor lasers and 
in designing even the complex structures of vertical cav- 
ity surface emitting lasers (VCSEL). One phenomenon, 
however, connected with VCSELs points to a deficiency 
of the SBE: these laser structures are known for their 
polarization instability, i.e., the uncontrolled switching 
of the laser modes between the two possible transverse 
polarizations. 0, Q, In addition, the investigation of 
semiconductor quantum structures as model systems for 
coupled Rabi oscillations with electrons, heavy- and 
light-hole states (each spin-degenerate) required to ex- 
tend the two-level SBE to six-level SBE and to take into 
account the polarization degree of freedom of the exciting 

light, a 

More recently, the carrier spin and its dynamics have 
gained much interest in the field of spintronics.|lO| 
Spin dynamics in semiconductors an d quantum 

structures, 0, 0] formulated so far in a more phc- 
nomenological way, is ruled by different mechanisms: 
one of which is related to the electron-hole exchange 
interaction. It becomes relevant if the semiconduc- 
tor system contains besides electrons also holes (e.g., 
due to doping or high excitation). This Bir-Aronov- 
Pikus (BAP) mechanism, originally considered for bulk 
semiconductors, has been described also for semi- 
conductor quantum structures. |l2l Il6| but never by a 
rigorous microscopic treatment of the spin dynamics. 
In this perspective, the SBE have been formulated for 
the six-level system, pj considering spin-splitting of the 
electronic energies due to spin-orbit coupling caused by 
bulk inversion (BIA) |17| or structure inversion asymme- 
try (SIA).[3 These extended SBE were designed only 
within the Hartree-Fock truncation leading to the coher- 
ent regime, thus neglecting scattering processes, respon- 
sible for relaxation. Recently, we have used the extended 



optical Bloch equations (SBE without carrier-carrier in- 
teraction |19|) to provide a microscopic approach to the 
longitudinal (2~i) and transverse ( Ta) relaxation times 
due to electron-phonon interaction. |20j In this approach 
we have considered scattering between electrons and 
phonons in second order Born approximation to provide 
a microscopic formulation for the D'yakonov Perel' (DP) 
mechanism of spin-relaxation. [2lJ The analogous concept 
is applied here to the electron-hole exchange interaction 
and yields the microscopic formulation of spin-relaxation 
due to the BAP mechanism. 

In the following, we concentrate on spin dynamics in a 
semiconductor quantum well (QW) under excitation with 
circularly polarized light leading to a nonequilibrium spin 
distribution due to optical selection rules. Let the system 
be described by the Hamiltonian 



(1) 



where 7io is the kinetic part including BIA and SIA spin- 
orbit coupling, 7iii g ht the interaction with the exciting 
light field, and 7i C oui the Coulomb interaction between 
the carriers. We adopt the notation of our previous works 
P, |2(j and use the basis in which the kinetic part of the 
Hamiltonian for the six-level system is diagonal 



k' m' 



Ho = J2 e <( k 0^(k')c<(k') 

e m ,(k')v m ,(k')vl,(k'). (2) 



Here, c mc (k) [v mv (k)] are fermion operators for elec- 
trons (light- and heavy-holes) with spin quantum num- 
bers m c = ±1/2 (m v = ±1/2 , ±3/2) defined with respect 
to the in-plane wave vector k. The time dependence of 
the operators is understood. The single particle energies 
e m '(k') \e m i (k')l describe subbands, which are spin-split 
due to spin-orbit interaction. In dipole approximation, 
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the interaction with the light field reads 

Wiight = - £ [E(<).d m , m ,(k')4,(k')^,(k') 



k' 



-E*(o-d; <m ,(k')« r< ,(k') Cl< (k') 



(3) 



where E(t) is the electric field vector and d^/ m i (k') is 
the dipole matrix element connecting valence and con- 
duction band states (for details see Ref. 0). 

The carrier-carrier interaction can be split up into four 
parts 



^coul — fLe 



IV i ivC i i_/X 

Hhh + n ch + n ch . 



(4) 



Here, TL ee (7ihh) describes the Coulomb interaction be- 
tween electrons (holes) in the conduction (valence) band. 
The remaining terms account for electron-hole interac- 
tion, the direct Coulomb term Ti^L and the exchange term 
Heh -E3 In tne f rame °f SBE.jjH especially for the 
coherent regime, carrier-carrier interaction has been con- 
sidered so far only with respect to renormalization of the 
single-particle energies and of the interaction with the 
light field, while the electron-hole exchange has been ig- 
nored. However, it is just H^ h which can cause spin-flips, 
thus contributing to the spin-dynamics due to the BAP 
mechanism and, hence, is of interest here. As derived in 
Ref. yj, the exchange term reads 



H oh = 2 E E m v m> c m> v (k, k', q) <£„ (~k + q) c r< (-k' + q) v} K (k) v mv (k') . (5) 

ra c m' c k k 



The detailed form of the interaction matrix element 
V* m m , m , (k, k', q) will not become important in the 
following. Rut we emphasize, that the structure of this 
matrix element makes simultaneous flips of electron and 
hole spins possible, which, in the electron system, finally 
contribute to spin- relaxation. 

While the dynamics of the whole system is contained 
in the equations of motion (EOM) of the 6x6 density 
matrix, we concentrate here on the dynamics of the elec- 



tron spins by looking at the EOM of the 2x2 density 
matrix for the electron subsystem 

g(m c fh c ) _ { £ra„m c (k) £?m c — m c (k) \ ^ 
y Q— m c m c (k) g„ m< ,_ me (k) J 

The single entries are expectation values of products of 
a creation and an annihilation operator g mc fn c (k) = 
(4, c (k) c^ c (k)} . Their EOM read 



i hd, 



6>m c m c (k) = fe mc (k) - e„ c (k) 



E E [■ 

k q m c 



V, 
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(k) + J2 [E(*) • d™ mv Pm c m v (k) - E* (t) • d™ c * mv Pt c mv (k) 

' rn> (~k + q, k, q) (c* (k) cw (-k + q) vt (-k + q) v M (k)) 



A „ mc ( k ' - k + q. q) ( c L (~ k + q) c ^ ( k ) vi, (k) Vrhv (-k + q)> , (7) 



where we have introduced the interband polarization 
P mcmv (k) = (4 lc (k) u^ v (k)).0 Due to the many-body 
contributions, the dynamics of g nlr r - nc (k) are ruled by 
four-point density matrices and, consequently, we run 
into a hierarchy problem^ which can be solved by an ap- 
propriate truncation. The Hartree-Fock (HF) truncation 
scheme Q, as used in Ref. Q, factorizes the expectation 
values of the four-operator terms into a product of two- 
operator terms under the condition that they are macro- 
scopic, namely, either electron (hole) densities or polar- 
izations. While closing the hierarchy and renormalizing 



the eigenenergies and the dipole interaction the HF trun- 
cation limits the EOM to the coherent regime, because 
no scattering processes are taken into account. In order 
to include these processes, which are essential for spin- 
relaxation and -dephasing, we go beyond the HF trunca- 
tion by considering the reduced four-operator terms 
defined as the difference between the expectation value of 
the untruncated four-operator term and its HF truncated 
product. For (c^, (-k + q) c„ c (k) yt , (k) u„„ (-k + q)) 
[see Eq. Q] it reads 
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6 ( c l>_ (- £ + q) c '»c (k) 4 1 (k) Vfh v (-k + q)) = (4, (-k + q) c rnc (k) 4, (k) tfe„ (-k + q)> 

-( C ^(k)c„- lc (k))( V t, (-k + q) UA „(-k + q))<S ki _ fc+q . (8) 



The scattering contributions are found by solving the 
EOM of the reduced four-operator terms which contain 
the complete information about the scattering in expec- 
tation values of four- and six-operator terms. In analogy 
to the case of electron-phonon scattering [2(| we trun- 
cate these terms by factorizing them into their macro- 
scopic parts and taking into account only those, which 
contribute in second order Born approximation. After 
integrating the arising equations and applying the adi- 
abatic and the Markov approximation |22(, we achieve 
a closed set of equations for the reduced four-operator 
terms, which can be solved and used in Eq. (for 
technical details see Ref. |2(|. Thus, the EOM for the 
diagonal entries of the 2x2 density matrix due to the 
electron-hole exchange scattering can be cast into the 



form 

d t Qm cmc Mh = ir*£ e (k)ft n .m (k) 

+C r X mc (k)[l-g mc „ lc (k)] , (9) 

with r™ tx c (k) [r™* mc (k)] accounting for the exchange 
scattering out of (into) the state with spin m c at wave 
vector k. The derivation of the scattering contribu- 
tions for the different four-operator terms in Eq. Q 
follows the same scheme. Thus, we present here only 
the results for the reduced four-operator term (c m i (— k + 

q) c TOc (k) v m , (k) Vm v (— k + q)). The corresponding out- 
scattering rate r™ tx (k) reads 



r™M(k) = J2 |Vm c ™;<m,(-k + q, k, q)| 2 5[e mc (k)-e^(k-q)- em ,(-k + q) + 6 r< (-k)] 

k q m' c 

X [1 - Qm' v m' v (~k)] [1 - g m ' a m' e (~k + q)] Qfh v m„ (k - q) Q mc rn c (k) , (10) 



with a similar expression for the in-scattering rate 
x m (k) . These expressions represent all contributions 
to electron-hole scattering by Coulomb exchange interac- 
tion in second order Born approximation. It is important 
to note that without a macroscopic occupation of hole 



states (by doping or optical excitation) this scattering 
rate vanishes: holes are required for the mutual spin flips 
of the BAP mechanism. 

The EOM for the off-diagonal entry of the density ma- 
trix the expressions can be written in the form 



<9 t £> mc -„ lc (k)|x = — 
* n 



m c — m 



(k) Qm c —m c (k) 



k q m c 



; (k + q) p mc _ mc (k + q) 



(11) 



As for the electron-phonon scattering [2(j, the first self- 
energy term in Eq. (|ll(l is proportional to the ab- 
solute squared value of the interaction matrix element 
"^m c ml m> m v ( — k + q, k, q) and can be split up into real 



and imaginary part connected by Kramers- Kronig trans- 
formation, where the imaginary part 
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3{S^ c _ mc (k)} = -J2 l V m c m;,^™„(- k + q> £ , q)| 2 <J[e m /(-k + q) - e_ mc (k) + e A „ (-k + q) -e A; (k)] 

kq raj. 

x | [1 - ftn; m; (-k)] [1 - Q m ' c m / (-k + q)] ^ (k - q) 

+ [1 - £>m„ m„ (k - q)] Qm' v rh' v (~k) £> m ^ m / (-k + q) } (12) 



accounts for dephasing due to scattering, while the real 
part contributes to the renormalization of the eigenener- 
gies. However, the real and imaginary part of the second 
term £^j c (k + q) in Eq. ljTl"|) are not connected by 
Kramers-Kronig theorem, because they are proportional 
to a product of two complex valued exchange interaction 



With the expressions given in Eqs. ©, (fTT)) . fT^ and 
(|13|1 we can define the T\ and T2 times of the Bloch 
equations as in Ref. l20t with the T\ time being ruled 
by the in- and out-scattering terms of Eq. © while the 
Ti time is given by the imaginary part of the self-energies 
of Eq. dj. The interpretation of these times with re- 
spect to spin-relaxation and -dephasing remains the same 
as in Ref. 20 apart from the fact that the spin-relaxation 
mechanism is different. 

In conclusion, we have derived microscopic expressions 
for the scattering rates due to electron-hole exchange in- 
teraction in a semiconductor QW in the frame of ex- 
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matrix elements. Nevertheless it is possible to sort out 
two parts, one proportional to principal values and one 
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der to derive the expressions corresponding to the ones 
given in Ref. l20t we will present here the part propor- 
tional to (5-functions denoted as f * c _ TOc (k + q) 



(13) 

I 

tended SBE. As it turns out, the expressions for these 
rates show the same qualitative structure as found for 
carrier-phonon scattering |20| . The particular scattering 
mechanism considered here is the one responsible for the 
BAP mechanism of spin relaxation. Thus, the presented 
results are a microscopic formulation of the BAP spin 
relaxation in the frame of the extended semiconductor 
Bloch equations. 
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